Ideas:

1 Transects

Plant Flowers Date lon lat ele Month Year julian
Glossoloma oblongicalyx 4 2015-10-19 -78.59093 0.130838 2270 October 2015 292
Gasteranthus quitensis 2 2016-10-17 -78.59770 0.120070 1940 October 2016 291
Kohleria affinis 1 2016-12-13 -78.59534 0.126746 2110 December 2016 348
Columnea ciliata 3 2014-02-27 -78.59934 0.116682 1960 February 2014 58
Columnea medicinalis 1 2014-04-23 -78.59372 0.128700 2130 April 2014 113
Drymonia teuscheri 3 2016-07-28 -78.59245 0.129393 2200 July 2016 210

2 Interactions

NEEDS TO BE FIXED.

3 Phylogeny

4 Traits

4.0.1 Total Flowers

4.1 Peak date

As range

4.2 Species elevation ranges

4.3 Flowering Data Matrix

5 Baseline model

Equal probability of flowering at anytime.

## sink("model/threshold_baseline.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
## 
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <- alpha[Plant[x]] 
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<-alpha[PredPlant[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
## 
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2548
##    Unobserved stochastic nodes: 3164
##    Total graph size: 15848
## 
## Initializing model

6 Get Chains

6.0.1 Evaluate convergence

6.0.2 Posterior estimates

7 Phylogeny

7.1 Attraction

## sink("model/threshold_attraction.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <-  alpha[Plant[x]] + e[Plant[x],Site[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Sites){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
## 
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=inverse(vCov*gamma)
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma = 1
##     
##     #Strength of covariance decay
##     lambda_cov = 3
##     omega  = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2548
##    Unobserved stochastic nodes: 3170
##    Total graph size: 19893
## 
## Initializing model

8 Get Chains

8.0.1 Evaluate convergence

8.0.2 Posterior estimates

8.1 Effect of phylogenetic attraction

8.2 Repulsion

## sink("model/threshold_repulsion.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <-  alpha[Plant[x]] + e[Plant[x],Site[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Sites){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
## 
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=vCov*gamma
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma = 5
##     
##     #Strength of covariance decay
##     lambda_cov = 2
##     omega = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2548
##    Unobserved stochastic nodes: 3170
##    Total graph size: 19893
## 
## Initializing model

9 Get Chains

9.0.1 Evaluate convergence

9.0.2 Posterior estimates

9.1 Effect of phylogenetic repulsion

10 Traits

10.1 Trait Attraction

## sink("model/threshold_attraction.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <-  alpha[Plant[x]] + e[Plant[x],Site[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Sites){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
## 
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=inverse(vCov*gamma)
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma = 1
##     
##     #Strength of covariance decay
##     lambda_cov = 3
##     omega  = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2548
##    Unobserved stochastic nodes: 3170
##    Total graph size: 20049
## 
## Initializing model

11 Get Chains

11.0.1 Evaluate convergence

11.0.2 Posterior estimates

11.1 Effect of trait attraction

11.2 Repulsion

## sink("model/threshold_repulsion.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <-  alpha[Plant[x]] + e[Plant[x],Site[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Sites){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
## 
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=vCov*gamma
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma = 5
##     
##     #Strength of covariance decay
##     lambda_cov = 2
##     omega = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2548
##    Unobserved stochastic nodes: 3170
##    Total graph size: 20049
## 
## Initializing model

12 Get Chains

12.0.1 Evaluate convergence

12.0.2 Posterior estimates

12.1 Effect of trait repulsion

13 Interaction

13.1 Attraction

## sink("model/threshold_attraction.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <-  alpha[Plant[x]] + e[Plant[x],Site[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Sites){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
## 
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=inverse(vCov*gamma)
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma = 1
##     
##     #Strength of covariance decay
##     lambda_cov = 3
##     omega  = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2548
##    Unobserved stochastic nodes: 3170
##    Total graph size: 20049
## 
## Initializing model

13.1.1 Evaluate convergence

13.1.2 Posterior estimates

13.2 Effect of interaction attraction

13.3 Repulsion

## sink("model/threshold_repulsion.jags")
## cat("
##     model {
##     
##     for (x in 1:Nobs){
##     
##     #Observation of a flowering plant
##     Y[x] ~ dbern(p[x])
##     logit(p[x]) <-  alpha[Plant[x]] + e[Plant[x],Site[x]]
##     
##     #Residuals
##     discrepancy[x] <- abs(Y[x] - p[x])
##     
##     #Assess Model Fit
##     Ynew[x] ~ dbern(p[x])
##     discrepancy.new[x]<-abs(Ynew[x] - p[x])
##     }
##     
##     
##     #Sum discrepancy
##     fit<-sum(discrepancy)/Nobs
##     fitnew<-sum(discrepancy.new)/Nobs
##     
##     #Prediction
##     
##     for(x in 1:Npreds){
##     #predict value
##     
##     #Observation - probability of flowering
##     prediction[x] ~ dbern(p_new[x])
##     logit(p_new[x])<- alpha[NewPlant[x]] + e[NewPlant[x],NewSite[x]]
##     
##     #predictive error
##     pred_error[x] <- abs(Ypred[x] - prediction[x])
##     }
##     
##     #Predictive Error
##     fitpred<-sum(pred_error)/Npreds
##     
##     #########################
##     #autocorrelation in error
##     #########################
##     
##     #For each of observation
##     for(y in 1:Sites){
##     e[1:Plants,y] ~ dmnorm(zeros,tauC[,])
##     }
## 
##     ##covariance among similiar species
##     for(i in 1:Plants){
##     for(j in 1:Plants){
##     C[i,j] = exp(-lambda_cov * D[i,j])
##     }
##     }
##     
##     ## Covert variance to precision for each parameter, allow omega to shrink to identity matrix
##     vCov = omega*C[,] + (1-omega) * I
##     tauC=vCov*gamma
##     
##     #Priors
##     
##     #Species level priors
##     
##     for (j in 1:Plants){
##     
##     #Intercept flowering probability
##     alpha[j] ~ dnorm(0,0.386)
##     
##     } 
##     #Autocorrelation priors
##     gamma = 5
##     
##     #Strength of covariance decay
##     lambda_cov = 2
##     omega = 1
##     }
##     ",fill=TRUE)
## 
## sink()
## Compiling model graph
##    Resolving undeclared variables
##    Allocating nodes
## Graph information:
##    Observed stochastic nodes: 2548
##    Unobserved stochastic nodes: 3170
##    Total graph size: 20049
## 
## Initializing model

14 Get Chains

14.0.1 Evaluate convergence

14.0.2 Posterior estimates

14.1 Effect of interaction repulsion

15 Model Comparison

15.0.1 E: The effect of autocorrelation

16 P(Flowering)

# One example

Glossoloma purpureum

17 By site

17.1 1300m - 1500m

17.2 1500-1700m

17.3 1700-1900

17.4 1900-2100m

17.5 2100-2300m

17.6 2300m - 2500m

17.7 Alpha: Species specific flowering rate

17.8 Omega: The magnitude of the effect of autocorrelation on mean flowering occurrence

17.9 Gamma: The variance of the effect of autocorrelation on mean flowering occurrence

17.10 Lambda: The decay in autocorrelation effect

17.11 Decay in autocorrelation effect

18 Model Fit

18.1 Bayesian pvalue

## # A tibble: 7 x 2
##   Model                       p
##   <chr>                   <dbl>
## 1 baseline                0.528
## 2 interaction_attraction  0.974
## 3 interaction_repulsion   0.978
## 4 phylogenetic_attraction 0.994
## 5 phylogenetic_repulsion  0.98 
## 6 trait_attraction        0.978
## 7 trait_repulsion         0.958

Without baseline 

## # A tibble: 6 x 2
##   Model                       p
##   <chr>                   <dbl>
## 1 interaction_attraction  0.974
## 2 interaction_repulsion   0.978
## 3 phylogenetic_attraction 0.994
## 4 phylogenetic_repulsion  0.98 
## 5 trait_attraction        0.978
## 6 trait_repulsion         0.958

18.2 Model Fit

Model mean lower upper
trait_attraction 0.2059569 0.1995477 0.2119726
phylogenetic_attraction 0.2046194 0.1980780 0.2109580
interaction_attraction 0.2043841 0.1977721 0.2109840
phylogenetic_repulsion 0.1973214 0.1911837 0.2035766
interaction_repulsion 0.1971715 0.1904616 0.2031607
trait_repulsion 0.1964914 0.1900542 0.2028746

18.2.1 Without baseline

18.3 By Species

18.3.1 Without baseline

18.3.2 Zoom in

19 Prediction

19.0.1 Tables

Model mean lower upper
baseline 0.2434784 0.2192691 0.2674419
phylogenetic_attraction 0.2126379 0.1910299 0.2358804
trait_attraction 0.2124917 0.1926080 0.2358804
interaction_attraction 0.2110764 0.1877076 0.2342193
interaction_repulsion 0.2098140 0.1876246 0.2358804
trait_repulsion 0.2097043 0.1877076 0.2342193
phylogenetic_repulsion 0.2084718 0.1843854 0.2342193

19.0.2 Flowering Rates among models

Plant Month Observed baseline interaction_attraction interaction_repulsion phylogenetic_attraction phylogenetic_repulsion trait_attraction trait_repulsion
Besleria solanoides Apr 13.3 17.0 15.7 16.9 16.7 16.8 16.1 16.3
Besleria solanoides Aug 5.9 16.6 16.5 16.5 16.4 16.5 15.6 17.7
Besleria solanoides Dec 12.5 17.8 16.6 17.6 15.9 18.3 16.1 17.9
Besleria solanoides Feb 0.0 16.6 16.0 16.6 15.5 16.0 15.8 16.0
Besleria solanoides Jan 0.0 17.3 16.8 18.1 17.6 17.1 17.9 16.3
Besleria solanoides Jul 5.6 16.8 14.9 15.3 15.0 15.8 14.6 14.8
Besleria solanoides Jun 5.0 16.1 15.2 15.8 15.5 15.4 15.5 16.6
Besleria solanoides Mar 0.0 17.4 16.8 17.4 17.5 17.1 17.4 17.3
Besleria solanoides May 0.0 16.1 16.5 17.4 17.2 16.5 16.4 17.1
Besleria solanoides Nov 33.3 16.4 16.1 17.5 17.3 15.6 15.9 17.6
Besleria solanoides Oct 55.6 16.7 16.7 17.8 17.9 16.9 17.4 16.9
Besleria solanoides Sep 50.0 15.5 16.3 15.8 15.5 16.4 15.5 15.7
Columnea ciliata Apr 46.7 15.8 17.8 18.8 17.8 19.3 17.4 18.9
Columnea ciliata Aug 0.0 16.6 13.5 12.7 14.2 13.0 13.9 12.4
Columnea ciliata Dec 0.0 15.4 16.2 15.3 16.7 15.3 14.6 15.2
Columnea ciliata Feb 43.8 16.2 15.5 16.5 16.2 16.2 15.8 15.6
Columnea ciliata Jan 25.0 15.9 19.0 19.2 18.0 19.6 18.3 19.3
Columnea ciliata Jul 5.6 16.1 16.4 17.4 16.8 17.1 16.7 16.4
Columnea ciliata Jun 10.0 15.7 15.8 17.0 15.9 15.7 16.2 15.5
Columnea ciliata Mar 66.7 16.1 17.9 19.4 17.8 19.9 18.3 20.2
Columnea ciliata May 5.9 16.6 15.9 17.0 15.5 16.9 16.8 16.6
Columnea ciliata Nov 0.0 17.1 16.6 16.4 15.8 16.6 16.0 15.5
Columnea ciliata Oct 0.0 16.4 15.7 14.8 15.1 15.3 15.2 14.5
Columnea ciliata Sep 0.0 15.7 15.6 16.3 16.2 16.6 16.1 15.8
Columnea kucyniakii Apr 6.7 8.8 8.6 8.9 7.9 8.9 8.4 8.2
Columnea kucyniakii Aug 0.0 8.8 8.3 9.3 8.4 9.4 8.8 9.2
Columnea kucyniakii Dec 12.5 8.5 9.0 9.5 9.8 8.9 9.3 9.5
Columnea kucyniakii Feb 43.8 8.1 8.2 9.6 8.4 8.5 8.4 9.1
Columnea kucyniakii Jan 25.0 8.9 9.5 9.3 8.4 10.2 9.2 9.8
Columnea kucyniakii Jul 0.0 8.7 7.9 8.8 7.6 8.3 8.2 8.4
Columnea kucyniakii Jun 0.0 9.2 8.1 8.7 8.9 8.5 8.1 8.8
Columnea kucyniakii Mar 16.7 8.8 8.4 10.2 9.0 9.8 9.0 8.9
Columnea kucyniakii May 0.0 8.4 8.7 9.6 9.3 9.6 10.1 8.8
Columnea kucyniakii Nov 11.1 8.4 8.5 8.1 7.8 8.5 7.6 9.6
Columnea kucyniakii Oct 0.0 8.4 8.6 8.6 8.3 9.3 7.9 9.1
Columnea kucyniakii Sep 0.0 8.8 8.2 9.6 8.2 8.9 9.1 8.9
Columnea medicinalis Apr 0.0 15.7 14.0 14.8 14.1 14.3 14.0 14.6
Columnea medicinalis Aug 23.5 14.1 16.8 16.0 15.6 16.6 16.2 15.6
Columnea medicinalis Dec 0.0 14.6 16.4 16.2 15.8 16.1 16.6 15.9
Columnea medicinalis Feb 12.5 14.0 14.2 14.5 14.1 14.1 14.1 14.2
Columnea medicinalis Jan 16.7 13.8 13.9 15.0 15.3 14.7 14.0 15.5
Columnea medicinalis Jul 27.8 14.3 13.2 13.0 13.2 13.0 12.4 13.2
Columnea medicinalis Jun 20.0 14.6 14.5 14.3 14.7 14.9 14.9 14.6
Columnea medicinalis Mar 8.3 14.3 14.0 15.2 14.8 15.4 13.4 15.0
Columnea medicinalis May 11.8 14.4 16.2 16.8 16.3 15.5 15.2 16.6
Columnea medicinalis Nov 11.1 14.1 15.5 14.6 14.6 15.5 15.1 14.7
Columnea medicinalis Oct 22.2 13.9 15.8 15.8 14.5 16.1 16.1 15.1
Columnea medicinalis Sep 5.0 13.8 14.3 14.6 14.5 14.1 14.0 14.5
Columnea picta Apr 13.3 14.0 14.4 16.0 14.8 15.8 14.4 15.6
Columnea picta Aug 5.9 14.3 12.2 11.3 12.1 11.5 11.8 11.4
Columnea picta Dec 12.5 14.3 14.7 14.4 15.0 13.6 13.9 13.9
Columnea picta Feb 6.2 14.2 14.0 15.6 14.5 14.8 14.6 13.9
Columnea picta Jan 16.7 14.5 14.6 14.8 15.1 14.7 14.7 15.4
Columnea picta Jul 16.7 14.0 14.0 16.1 15.2 16.0 14.9 15.2
Columnea picta Jun 25.0 14.5 14.0 14.8 13.8 14.4 13.9 14.3
Columnea picta Mar 0.0 14.8 14.5 14.9 14.2 15.4 14.0 14.6
Columnea picta May 35.3 14.1 17.1 17.9 18.0 17.6 17.4 17.4
Columnea picta Nov 0.0 13.9 10.6 10.6 10.5 10.3 9.1 10.4
Columnea picta Oct 0.0 13.7 10.4 9.1 9.5 9.5 9.6 9.6
Columnea picta Sep 20.0 14.0 15.5 16.2 15.9 16.8 15.7 16.5
Columnea strigosa Apr 0.0 12.2 11.5 11.5 11.1 10.8 11.6 11.4
Columnea strigosa Aug 11.8 12.3 14.6 15.0 14.2 15.6 15.3 15.4
Columnea strigosa Dec 0.0 11.5 12.2 11.2 11.6 11.3 12.7 11.3
Columnea strigosa Feb 12.5 12.6 12.4 13.9 12.7 13.2 13.1 13.7
Columnea strigosa Jan 8.3 11.6 9.2 10.1 9.5 9.5 9.8 9.5
Columnea strigosa Jul 27.8 12.7 12.8 14.3 13.9 13.8 13.3 14.7
Columnea strigosa Jun 40.0 12.0 14.0 14.0 12.7 13.1 13.9 13.8
Columnea strigosa Mar 0.0 12.0 8.9 10.1 9.8 9.3 10.0 9.2
Columnea strigosa May 23.5 12.1 10.4 9.9 10.1 9.6 10.2 9.7
Columnea strigosa Nov 0.0 12.3 13.9 14.9 13.4 15.3 14.2 14.3
Columnea strigosa Oct 0.0 12.6 14.5 14.2 13.2 14.6 14.3 15.0
Columnea strigosa Sep 0.0 11.7 12.4 12.0 11.9 11.9 12.4 12.7
Drymonia collegarum Apr 13.3 12.8 13.2 12.9 13.7 13.1 12.9 13.5
Drymonia collegarum Aug 17.6 13.2 14.9 14.9 16.1 14.2 15.6 13.8
Drymonia collegarum Dec 25.0 14.2 15.1 15.0 16.0 16.2 15.6 15.0
Drymonia collegarum Feb 6.2 13.3 13.1 12.8 12.5 12.4 13.2 12.4
Drymonia collegarum Jan 16.7 13.1 15.2 15.2 13.9 15.0 13.4 15.0
Drymonia collegarum Jul 11.1 13.0 12.0 11.5 12.5 10.9 12.0 11.5
Drymonia collegarum Jun 15.0 13.6 13.2 13.9 13.3 13.3 13.2 12.8
Drymonia collegarum Mar 8.3 13.2 15.8 15.2 14.1 15.7 13.2 15.6
Drymonia collegarum May 23.5 13.5 15.5 14.7 13.8 15.2 13.9 14.8
Drymonia collegarum Nov 11.1 14.3 13.7 14.4 14.8 13.4 14.3 13.5
Drymonia collegarum Oct 5.6 13.5 13.2 14.7 14.1 13.9 14.4 14.0
Drymonia collegarum Sep 10.0 13.9 13.3 12.9 13.0 13.1 12.4 13.0
Drymonia tenuis Apr 13.3 13.9 14.7 15.1 14.6 15.5 14.3 15.4
Drymonia tenuis Aug 11.8 14.1 11.0 13.7 13.3 11.8 13.6 11.6
Drymonia tenuis Dec 0.0 14.1 13.0 13.3 13.1 13.1 13.9 12.8
Drymonia tenuis Feb 18.8 14.1 14.8 14.0 14.6 14.9 13.2 14.5
Drymonia tenuis Jan 16.7 14.8 15.8 17.2 17.0 17.0 15.8 17.2
Drymonia tenuis Jul 16.7 13.9 13.8 13.8 13.6 13.9 12.7 14.7
Drymonia tenuis Jun 20.0 14.6 13.8 13.0 13.8 13.7 13.2 13.8
Drymonia tenuis Mar 25.0 14.2 17.2 17.0 17.1 17.0 15.2 17.5
Drymonia tenuis May 23.5 13.4 14.6 13.1 13.9 12.5 13.5 13.1
Drymonia tenuis Nov 0.0 14.8 13.6 16.1 15.5 14.8 15.6 16.0
Drymonia tenuis Oct 5.6 14.3 13.9 16.5 15.3 15.2 16.0 14.5
Drymonia tenuis Sep 5.0 13.5 13.3 13.3 13.4 13.7 13.3 13.5
Drymonia teuscheri Apr 26.7 16.4 15.9 15.4 16.3 16.1 16.4 16.6
Drymonia teuscheri Aug 23.5 16.5 19.3 19.7 19.4 20.1 18.8 20.7
Drymonia teuscheri Dec 0.0 16.4 17.3 18.1 17.3 17.4 18.2 18.0
Drymonia teuscheri Feb 6.2 16.7 16.5 16.8 16.8 16.9 17.2 17.9
Drymonia teuscheri Jan 16.7 16.7 12.7 14.4 14.1 13.2 15.7 13.3
Drymonia teuscheri Jul 22.2 17.5 17.3 16.3 15.9 17.8 16.7 17.6
Drymonia teuscheri Jun 15.0 17.0 17.4 17.4 17.5 18.2 17.1 17.5
Drymonia teuscheri Mar 16.7 16.7 13.2 14.4 15.1 13.2 14.1 13.2
Drymonia teuscheri May 29.4 17.0 14.7 16.1 16.2 15.1 16.0 14.9
Drymonia teuscheri Nov 11.1 16.3 18.8 18.8 18.1 18.9 18.1 20.0
Drymonia teuscheri Oct 11.1 16.6 18.8 19.1 18.8 19.8 18.6 19.2
Drymonia teuscheri Sep 10.0 16.6 16.8 16.4 16.6 17.1 16.2 16.8
Gasteranthus lateralis Apr 13.3 14.0 13.3 14.4 13.5 13.7 14.8 13.8
Gasteranthus lateralis Aug 17.6 14.0 16.0 16.0 15.6 16.4 15.1 16.8
Gasteranthus lateralis Dec 25.0 14.6 15.2 15.0 14.5 14.1 15.3 15.8
Gasteranthus lateralis Feb 0.0 14.5 14.2 14.2 14.1 14.3 13.9 14.1
Gasteranthus lateralis Jan 33.3 14.8 13.0 12.8 12.9 12.9 14.1 12.2
Gasteranthus lateralis Jul 0.0 13.6 13.0 14.6 13.5 13.9 13.1 13.8
Gasteranthus lateralis Jun 0.0 14.3 14.2 15.2 13.9 14.8 13.9 14.6
Gasteranthus lateralis Mar 8.3 13.7 13.8 13.3 12.8 12.7 13.7 11.8
Gasteranthus lateralis May 5.9 13.3 13.6 14.1 14.4 14.0 13.7 13.0
Gasteranthus lateralis Nov 44.4 14.5 14.8 14.8 15.0 16.1 14.6 15.7
Gasteranthus lateralis Oct 33.3 14.0 14.6 15.8 14.8 15.4 15.7 15.4
Gasteranthus lateralis Sep 10.0 14.4 13.4 14.4 14.3 14.0 14.0 14.3
Gasteranthus quitensis Apr 0.0 11.7 12.1 13.0 12.3 12.4 11.6 12.2
Gasteranthus quitensis Aug 17.6 12.3 11.2 10.0 10.3 10.0 10.4 11.0
Gasteranthus quitensis Dec 12.5 11.9 11.4 10.9 10.2 11.0 11.5 11.5
Gasteranthus quitensis Feb 18.8 12.0 11.4 13.8 12.3 12.4 12.7 12.5
Gasteranthus quitensis Jan 33.3 12.3 12.8 14.4 14.2 13.1 14.1 14.1
Gasteranthus quitensis Jul 5.6 11.7 11.8 12.7 11.8 12.0 11.7 12.9
Gasteranthus quitensis Jun 0.0 11.4 11.1 11.8 11.2 11.3 10.8 11.5
Gasteranthus quitensis Mar 0.0 11.5 13.0 14.8 14.0 13.3 13.8 14.6
Gasteranthus quitensis May 0.0 12.2 11.6 11.5 11.6 11.0 11.3 12.4
Gasteranthus quitensis Nov 11.1 12.5 12.2 11.5 12.5 10.8 11.6 11.1
Gasteranthus quitensis Oct 22.2 12.1 12.9 12.3 12.8 12.0 12.5 12.1
Gasteranthus quitensis Sep 20.0 12.2 11.9 12.2 11.5 11.6 11.9 13.1
Glossoloma oblongicalyx Apr 0.0 18.7 17.3 17.0 17.5 17.1 17.5 16.8
Glossoloma oblongicalyx Aug 35.3 18.6 21.9 22.6 21.4 22.7 21.2 23.5
Glossoloma oblongicalyx Dec 0.0 18.0 18.2 18.0 17.8 18.3 18.9 18.4
Glossoloma oblongicalyx Feb 0.0 18.1 18.6 20.2 19.5 19.4 18.9 19.4
Glossoloma oblongicalyx Jan 0.0 17.4 14.4 13.3 14.6 14.0 14.4 13.4
Glossoloma oblongicalyx Jul 38.9 18.6 18.9 21.0 19.6 19.7 20.1 20.7
Glossoloma oblongicalyx Jun 10.0 18.4 19.8 20.2 20.3 20.2 19.6 20.1
Glossoloma oblongicalyx Mar 0.0 19.0 14.6 14.0 13.7 14.7 14.7 13.6
Glossoloma oblongicalyx May 11.8 18.4 15.4 14.1 15.7 14.6 16.5 14.8
Glossoloma oblongicalyx Nov 22.2 18.4 20.5 21.0 20.2 22.9 19.7 23.2
Glossoloma oblongicalyx Oct 38.9 18.5 21.1 21.4 20.5 21.9 20.4 22.9
Glossoloma oblongicalyx Sep 35.0 18.7 17.8 18.6 18.1 17.6 17.9 18.0
Glossoloma purpureum Apr 6.7 13.4 14.8 15.7 15.4 14.4 15.5 15.4
Glossoloma purpureum Aug 17.6 13.9 10.3 10.6 10.8 10.5 12.1 10.0
Glossoloma purpureum Dec 12.5 13.1 13.2 13.4 12.6 13.9 13.4 12.5
Glossoloma purpureum Feb 12.5 13.2 13.7 13.7 13.8 13.5 14.0 13.5
Glossoloma purpureum Jan 16.7 14.3 14.1 14.1 13.2 15.2 13.9 13.9
Glossoloma purpureum Jul 22.2 13.7 14.8 15.4 15.2 14.9 15.9 14.2
Glossoloma purpureum Jun 10.0 12.9 13.8 14.0 14.0 13.8 14.3 14.0
Glossoloma purpureum Mar 0.0 13.4 13.9 13.9 14.0 14.8 13.3 14.3
Glossoloma purpureum May 5.9 12.7 16.3 18.7 16.8 17.8 16.9 17.9
Glossoloma purpureum Nov 11.1 14.6 10.7 8.9 10.0 8.8 10.8 8.6
Glossoloma purpureum Oct 11.1 13.2 9.9 8.5 9.3 8.7 10.4 8.3
Glossoloma purpureum Sep 25.0 13.9 15.2 15.6 15.2 15.2 15.0 15.3
Kohleria affinis Apr 46.7 15.5 15.9 15.1 15.1 16.4 15.6 15.8
Kohleria affinis Aug 5.9 15.0 16.2 16.6 17.0 16.2 17.1 15.7
Kohleria affinis Dec 12.5 15.8 17.6 17.3 17.6 16.4 17.8 16.9
Kohleria affinis Feb 25.0 15.4 15.0 14.4 15.6 14.8 15.1 15.4
Kohleria affinis Jan 8.3 14.9 16.6 16.7 15.2 15.7 16.6 16.1
Kohleria affinis Jul 0.0 15.4 14.4 14.0 13.8 15.5 14.3 14.2
Kohleria affinis Jun 15.0 15.3 14.7 14.6 15.8 15.8 15.2 15.6
Kohleria affinis Mar 41.7 16.1 16.1 16.0 16.5 15.6 15.8 15.6
Kohleria affinis May 29.4 15.7 15.5 16.3 17.2 16.8 16.5 15.4
Kohleria affinis Nov 0.0 15.4 16.8 15.6 17.5 15.5 17.2 17.2
Kohleria affinis Oct 0.0 14.9 17.1 16.6 17.3 16.2 16.8 16.4
Kohleria affinis Sep 5.0 15.1 15.1 15.5 15.4 14.9 15.6 15.8

20 Flowering correlations in predictions

Grab an example Drymonia teuscheri v Drymonia collegarum

D["Drymonia teuscheri","Drymonia collegarum"]
## [1] 0.02974
Dtraits["Drymonia teuscheri","Drymonia collegarum"]
## [1] 0.04191869
Dint["Drymonia teuscheri","Drymonia collegarum"]
## [1] 0.4755805

Logit E

inv.logit E 

## # A tibble: 6 x 4
## # Groups:   Model [6]
##   Model                   Var1               Var2            Correlation_E
##   <chr>                   <fct>              <fct>                   <dbl>
## 1 interaction_attraction  Drymonia teuscheri Drymonia colle…       -0.492 
## 2 interaction_repulsion   Drymonia teuscheri Drymonia colle…        0.0870
## 3 phylogenetic_attraction Drymonia teuscheri Drymonia colle…        0.784 
## 4 phylogenetic_repulsion  Drymonia teuscheri Drymonia colle…       -0.864 
## 5 trait_attraction        Drymonia teuscheri Drymonia colle…        0.698 
## 6 trait_repulsion         Drymonia teuscheri Drymonia colle…       -0.867
## # A tibble: 7 x 4
## # Groups:   Model [7]
##   Model                   Var1               Var2            Correlation_P
##   <chr>                   <fct>              <fct>                   <dbl>
## 1 baseline                Drymonia teuscheri Drymonia colle…        NA    
## 2 interaction_attraction  Drymonia teuscheri Drymonia colle…        -0.394
## 3 interaction_repulsion   Drymonia teuscheri Drymonia colle…        -0.160
## 4 phylogenetic_attraction Drymonia teuscheri Drymonia colle…         0.495
## 5 phylogenetic_repulsion  Drymonia teuscheri Drymonia colle…        -0.469
## 6 trait_attraction        Drymonia teuscheri Drymonia colle…         0.355
## 7 trait_repulsion         Drymonia teuscheri Drymonia colle…        -0.471
## # A tibble: 7 x 4
## # Groups:   Model [7]
##   Model                   Var1               Var2         Correlation_Ynew
##   <chr>                   <fct>              <fct>                   <dbl>
## 1 baseline                Drymonia teuscheri Drymonia co…           0.0889
## 2 interaction_attraction  Drymonia teuscheri Drymonia co…          -0.394 
## 3 interaction_repulsion   Drymonia teuscheri Drymonia co…          -0.162 
## 4 phylogenetic_attraction Drymonia teuscheri Drymonia co…           0.490 
## 5 phylogenetic_repulsion  Drymonia teuscheri Drymonia co…          -0.468 
## 6 trait_attraction        Drymonia teuscheri Drymonia co…           0.355 
## 7 trait_repulsion         Drymonia teuscheri Drymonia co…          -0.471